Brans-Dicke theory of gravity with torsion: A possible solution of the cj-problem 
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We study the Brans-Dicke theory of gravity in Riemann-Cartan space-times, and obtain general torsion solu- 
tions, which are completely determined by Brans-Dicke scalar field $, in the false vacuum energy dominated 
epoch. The substitution of the torsion solutions back to our action gives the original Brans-Dicke action with 
<I>-dependent Brans-Dicke parameter a;($). The evolution of during the inflation is studied and it is found 
that uj approaches to infinity at the end of inflation. This may solve the cj-problem in the extended inflation 
model. 

PACS numbers: 98.80.Cq, 98.80.Jk 



I. INTRODUCTION 

The discovery of spin- 1/2 fermions indicates that matter 
fields have microscopic structure, which can be character- 
ized by the spin angular momentum 3-forms Sab- Since gen- 
eral relativity (GR) is established (by hypothesis) in pseudo- 
Riemannian (i.e. torsion free) space-times, the gravitating 
sources are solely described by the stress-energy tensor, and 
intrinsic spin does not play any role in the conservation laws 
of angular momentum. Hence, GR lacks of a description of 
spin-orbital coupling. This unsatisfactory situation actually 
happens in any theory of gravitation established in pseudo- 
Riemannian space-time, and may be resolved when we extend 
relativistic theories of gravity to Riemann-Cartan (RC) space- 
time. Several well-kown theories of gravitation, e.g. Poincare 
gauge theory of gravity (PGT) (see a review article 1 1 1), are 
built in RC space-times. Moreover, it was discovered that a 
Lorentz gauge-covariant form of the Brans-Dicke (BD) the- 
ory of gravity yields torsion fields determined by the gradient 
of the BD scalar field $ EH. 

In RC space-times the concept of the metric g and metric- 
compatible connection V are fundamentally independent, so 
the associated independent variables of gravitational fields are 
orthonormal co-frames {e°} and connection 1-forms {i^°'b] 
[f4|. The intrinsic spin quantities Sab become the sources 
of torsion fields, and the Bianchi identities yield the conser- 
vation law of orbital angular momentum and intrinsic spin. 
Since torsion and the intrinsic spin have direct interactions, 
spin-polarized bodies can be used to detect torsion directly 
in the laboratory (see a review article |5|). Up to now, there 
is no experimental evidence showing the existence of torsion 
field, so the constraints on torsion-spin coupling turn out to 
be extremely small |5, 6|. The result is not surprising since 
it is very difficult to produce a significant magnitude of the 
intrinsic spin in the laboratory. However, we expect that in- 
trinsic spin should generate observable torsion effects in the 
early Universe. Refs |7 , 8 1 discovered that the spin-0+ torsion 
(sometimes called scalar torsion) mode in PGT |9| can nat- 
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urally explain the current acceleration of our Universe with- 
out introducing the dark energy. Moreover, we found that the 
quadratic curvature terms in RC space-times can generate a 
power-law inflation without introducing inflaton fields |10|. 
However, this inflation model is not satisfactory since it re- 
quires some fine-tuning in the parameters. 

The old inflation model was originally proposed to solve 
horizon and flatness problems by considering the Universe to 
undergo a first-order phase transition 1 1 1 1. Cosmic inflation is 
driven by the false vacuum energy of an inflaton field a and is 
supposed to end by bubble nucleations. It was soon realized 
that inflation will never come to an end because of the small- 
ness of the (dimensionless) bubble nucleation rate e = X/H^, 
where A is the number of bubbles nucleated per unit time per 
unit 3-voIume and H — a/a is the Hubble parameter during 
inflation lfT2ll . If we assume that bubble nucleation is domi- 
nated by quantum-mechanical tunneling, A can be expressed 
as A = Ae^^'^ in the semiclassical limit lfT3l - [T6l . The pref- 
actor A is equal to times terms expected to be of order 
unity and Se is the Euclidean action of the bounce solution, 
where Tc denotes the critical temperature of the phase transi- 
tion. Since the scale factor a{t) is exponentially expanding, 
a oc e^*, and if is a constant during inflation, we obtain that 
e is a constant. In lfT2]| . Guth and Weinberg proved that if 
e > e^r ~ 0.24, the system of bubbles will percolate at some 
finite time. However, a direct calculation of e from effective 
potentials of some specific models shows that the value of e 
is quite small (10^^"° is even plausible) [12|. It turns out that 
inflation never comes to the end in the old inflation model, and 
this is called the "graceful exit" problem. 

La and Steinhardt |17 | discovered that the "graceful exit" 
problem can naturally be resolved in the BD theory of gravity 
IfTSll . and they called this model the extended inflation. The 
main feature of the extended inflation is that a{t) has a power- 
law solution a oc t"+3 instead of exponential expansion, 
where w denotes the dimensionless BD parameter, so that 
e{t) cx is now time-dependent and monotonically increases 
with respect to time. It means that e can be very small in 
the beginning of inflation and then grows to the critical value 
Ecr, where the system of bubbles will percolate. Although 
the Universe can exit from false-vacuum energy domination 
in the extended inflationary scenario, it was soon realized that 
in order to satisfy the nearly isotropic spectrum of the cosmic 
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microwave background radiation (CMB), the constraint on a 
bubble-size distribution requires a; < 25 lfT9H2Tll . However, 
the current solar system observation of Cassini spacecraft re- 
quires that Lo must exceed 40000 1221 1231 . Apparently, the 
constraint of the bubble-size distribution is conflict with the 
solar system observations. This is called the w-problem. 

There are several phenomenological approaches to solve 
the cij-problem by assuming to be a function of $ Ii25 j or 
adding a potential V{^) in the BD action |l20l. In |I23, 
is taken to be = wq + Wm(T^)™, which increases 

monotonically during the evolution of $. Here, Mpi denotes 
the Planck mass and uj,n must be assigned a huge value to 
satisfy the solar system observations. When $ approaches to 
A/p; in the post-inflationary state, cj„i dominates, and the so- 
lution of a{t) gives a{t) cx ex^iatf), where a is a positive 
constant and f = t^™., . This solution is called the inter- 
mediate inflation. A further investigation of the bubble-size 
distribution in this intermediate inflation gives a constraint on 
cjQ and 771 |26|. Roughly speaking, it requires that the transi- 
tion between cjq and must be rapid, which corresponds to 
requiring a large value of m. 

The Einstein-Cartan theory is a natural generalization of 
GR to RC space-times, and its torsion fields are completely 
determined by the distribution of the intrinsic spin Sab- If the 
magnitude of Sab is too small to observe, there is no difference 
between the Einstein-Cartan theory and GR. However, this is 
not true in the BD theory with torsion, since torsion fields will 
be produced not only by the intrinsic spin but also by the gra- 
dient of BD scalar field $ |2J . Here, the torsion field generated 
by <i> is called BD torsion field. Hence both Sab and $ become 
the sources of torsion fields. In this paper, we find that the BD 
torsion field will contribute to cj and obtain the explicit form 
of cj($). Moreover, we show that the w-problem can naturally 
be resolved in the BD theory with torsion, instead of using 
phenomenological approaches. lo{^) only contains one phys- 
ical parameter 02, and the present value of $, which is equal 

Hence, w(<i>) 



to Mil 



Mpi (see Sec. 
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will yield 02 

actually does not contain any free parameter To understand 
whether this extended inflation model with torsion can satisfy 
the CMB anisotropic observation, it requires a further investi- 
gation on the cosmological perturbation in RC space-times. 

This paper is outlined as below. In Sec. [II] we generalize 
the BD theory to RC space-times. The main difference be- 
tween our work and Dereli & Tucker's work |2] is that our 
BD action with torsion includes three irreducible pieces of 
quadratic torsion, which were not considered in |2|. These 
quadratic torsion terms may be associated to kinematic en- 
ergy of orthonormal co-frames e°. Adding these terms does 
not spoil the field equations as the second-order differential 
equations. Without introducing any symmetry of space-time, 
we obtain a general torsion solution completely determined 
by $, where the Lagrangian of matter is assumed to be the 
potential U{a) of the inflaton field. In Sec. II A we substitute 
the torsion solution back to our original action, and obtain an 
effective action, which is equivalent to the original BD theory 
except that now is a function of $ instead of the dimen- 
sionless parameter. In Sec. Ill we study field equations in the 
homogeneous and isotropic Universe, and obtain analytic and 



numerical solutions of a{t) and during the inflation. Sec. 
JV] gives a discussion and conclusion. In this paper, we use 



the units c = ^ = 1 and SnG = Mp^ 



II. BRANS-DICKE THEORY OF GRAVITY WITH 
TORSION 

Mach's principle is a fundamental principle to explain the 
origin of inertia |28|. In attempting to incorporate Mach's 
principle, the BD theory introduces an inertial field $ which 
plays the role of the gravitational constant G and is deter- 
mined by the matter field distributions. So the gravitational 
fields are described by the metric g and the BD scalar field $, 
which has the dimension [$] = [Af]^. The BD theory starts 
from the following action: 

^BD($,e'^;*) = / I i?afc A *(e" A e') 



Wo 

2$ 



A^d^ + LmI^), (1) 



where Rab are Riemann curvature 2-forms, ujq is the BD di- 
mensionless parameter, and * is the Hodge map associated to 
4-dimensional metric g. denotes the Lagrangian 4-form 
of matter fields 'J, where the minimal coupling of gravita- 
tional fields is assumed, so there is no direct interaction be- 
tween vj/ and $. An important feature of BD theory is that 
when LjQ approaches to infinity, the field equations of $ yields 
that $ becomes a constant $0. Hence, the BD theory will 
recover to GR in the limit of 00. 

The most natural generalization of the BD theory to RC 
space-times is to start from the following action: 

3 

En n ujr, 

a„ T A * T a - ^ rf* A *d$ + Lm{^), (2) 



where Rab are curvature 2-forms defined by Rab = dujab + 

n 

(jjac A uj'^b, and T " denotes three irreducible pieces of torsion 

2-forms T"^ defined by T" = de" + uj°'c A e^. One may notice 

that the orthonormal co-frames e° and the connection 1 -forms 

ujab become independent variables. The trace vector torsion 
2 3 
T " (scalar torsion) and the axial torsion T " (pseudo-scalar 

torsion) can be expressed as 



2 



^fenAe^ 



T^-^li^iepATP), (3) 



where ia denotes the interior derivative, and the tensor part of 

1 1 2 3 

torsion T " is defined hy T = T"-- T T °-. Clearly, 
the dimensions of three parameters a„ are the same and equal 
to [an] = [Af]2. In Eq. we also assume the minimal 
coupling of the matter fields and the gravitational fields, so 
$ does not appear in the Lagrangian Lm- Generally speak- 
ing, the difference between Lm and Lm occurs when has a 
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direct interaction with LUat in the Lagrangian 4-forms, and ac- 
cording to the standard model of particle physics, only spin- 
1/2 fermions will directly couple to the connection 1 -forms 
in the action. This is the reason why fermions become the 
sources of torsion in RC space-times. In Eq. (|2]l, the first 
term on the right-hand side shows that $ couples to the full 
scalar curvature instead of the Riemannian scalar curvature, so 
one can expect that $ will generate torsion through these cou- 
pling. Furthermore, the field equations remain second-order 
differential equations when adding three irreducible quadratic 
torsion terms. In this paper, we will concentrate on the evolu- 
tion of $ and the gravitational fields at the inflationary epoch, 
when is dominated by the false vacuum energy, so we put 
Lm = —U{(j) * 1, where the potential U (cr) of the inflaton a 
is constant during inflation. 



Since the field equation 



SSe 



yields the algebraic 



equations for torsion 2-forms T°, we first solve these equa- 

2 

tions and then obtain general solutions T ° = e° A 2{<i>^a2) 

1 3 

with T " =T " = 0. One should notice that, to obtain 

1 3 

T " =T ° = 0, we have excluded degenerate situations where 
$ = ai or $ = 2a3. It is clear that the coupling of $ and the 

2 

scalar curvature will only produce scalar torsion T When 
02 = 0, we return to the result in [ Jl- The substitution of the 

n 

torsion solutions T " back to the remaining field equations, 
i.e. = and = 0, yields 



(4) 



+2D * Tc} = -^{icd<^ *d^ + d<^>Aic*d<^) + U* e^, 

$ . 1 

= i?afcA*e"''H d$Ad*$, (5) 

2(jJo 2$ 



where the BD torsion is 



rj-ia 



A 
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e" A 



2($-a2)' 



(6) 



A e" A A 



A Cd and D denotes 

1 3 

covariant exterior derivative II29II . Since T ° and T " vanish, 
we expect that Eq. Q does not contain the parameters oi 
and 03. It is clear that Eqs. Q and (j5]l are the second-order 
differential equations for e"^ and $. In Subsection II A we 



will show that the effects of the BD torsion can actually be 
combined with ujq to form an effective BD "parameter" 
which is now a function of $. More specifically, Eqs. Q 
and (|5]l are equivalent to original BD equations in Riemannian 
space-time with the effective BD parameter w(<i>). 



ujab associated with the Levi-Civita connection and the con- 
torsion 1 -forms Kab, which are defined by 1,291 

LJab = ^{ePiaibdep + ibdCa-iadeb), 
Kab = ^{-e^'iahdTp - itdTa + iadTb). (7) 
Substituting ujab ~ '^ab + Kai, into the definition of Rab yields 



Rah — Rn 



DKa 



Kar A 



(8) 



where t) is the covariant exterior derivative associated to iJjab- 
Using Eq. (|6]l, we obtain 

1 



Ka 



2($-a2) 



(ea ibd<^ - ebiad^). 



(9) 



The substitution of Eqs. (|8]l and (|9]l into Eq. Q gives the 
effective BD action: 

5Bi,(e",$;a) = I -Rab A ^e^' ^ ^d^ /\ ^d^ 

+C/(cr) * 1 + d6, (10) 
where the effective BD parameter is 
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2(02 - 



(11) 



and dB denotes the boundary term. If 02 — 0, we obtain 
uj = — |, which agrees with the result in [|2j|. It is not 
difficult to verify that the field equations obtained by varying 
Sbd with respect to e° and $ are equivalent to Eqs. ^ and 
(|5]l with the substitution of Eqs. ([8])-(^. 

Before we present a detail study of the evolution of 
and the scale factor a{t) during the inflation, we can first ex- 
amine the behavior of a;(<I>). Eq. (Ill indicates that when 
$ <C 02, we have w « wq- Moreover, when $ approaches to 
02, a;($) will then approach to infinity. It is clear that 
monotonically increases with respect to the growth of $. In 
order to satisfy the constraint uj > 40000 of the current so- 
lar system observation ll22ll . it requires that $ should be very 



close to 02 at present time. In Sec. Ill we apply analytic and 



numerical approaches to study the evolution of during 
the inflation, which gives that $ is asymptotically approach to 
02 in the post-inflationary stage. Furthermore, the equations 
of motion indicate that $ continuously approaches to 02 in the 
radiation and matter domination epochs, so this result can be 
used to explain why oj{^) is so large at the present time. 



III. EQUATIONS OF MOTION IN ROBERTSON-WALKER 
SPACE-TIMES 



A. The effective action 

In order to compare Eqs. (j4]) and (j5]l to the original BD 
equations, we should decompose the curvature 2-forms Rab 
into Riemannian curvature 2-forms Rab and torsion parts. The 
first step is to decompose LUab into the connection 1 -forms 



Although our Universe is apparently inhomogenous and 
anisotropic in small scales (e.g. galactic scale), the astro- 
physical observations strongly support the homogeneity and 
isotropy of our observable Universe in the cosmological scale. 
It allows us to assume that our observable Universe exists 3- 
dimensional space-like hypersurfaces, which are maximally 
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symmetric 3-spaces f3Q\. The assumption of homogeneity 
and isotropy in RC space-times gives 



e° = dt. 



T° = 0, T" 



a{t) 



il^lkr^)-^ ^^^^ 
f{t) e" A e" + h{t) * (e° A e"),(13) 



and $ = ^{t), where k — { — 1,0, 1} denotes the constant 
curvature of 3-dimensional spaces and r = y/x^Xa- It is con- 
venient to introduce a dimensionless scalar field x defined by 
X = ^ . The substitution of Eq. ( |l3| into Eq. (j6j yields 

X 



fit) 



2(X-1)' 



(14) 



and h{t) = 0. Moreover, the substitution of Eqs. (12i-(14i 
into Eqs. (|4])-(j5]l yields 

4x(l-x) 6 U 
3 3x' 



X 



Wo 



2(1 -X) 



(X + SHx) = - 



4(1 -X)^ 



/if 

:^^,(15) 
3xa2 

^(16) 



02 



where = and [/ = Afi ||3T|. Mp denotes the false 



vacuum energy, which may be around the GUT energy scale 
10^'* Gev. We should now try to determine the energy scale of 
a2. Since x(^) at present time tp is extremely close to 1, i.e. 
<I>(tp) w a2, and $(ip) should be normalized to (SttG)^^, 
we obtain that 02 ~ Mp^ . 

Eqs. ( 15 1-( 16 1 describe the evolution of a{t) and x(i) dur- 
ing the inflation. We first observe that Eq. ( [T6| has a very 
interesting feature. On the right-hand side of Eq. (I61, the 
first term is definitely negative and is proportional to x^, so 
one may identify it as frictional force. Moreover, the second 
term is a definitely positive constant, so it can be considered 
as a constant external force supplying x with kinetic energy. 
If we consider x ^ 1 at the beginning of inflation, the first 
term can actually be neglected so the false vacuum energy will 
drive x to have positive velocity and acceleration. It means 
that x{t) grows with positive acceleration. However, when x 
is approaching to 1, the frictional force cannot be neglected 
anymore. So one can expect that x{t) will evolve from accel- 
erating phase to decelerating phase. In Subsection III A we 



16 1 in the early 



obtain analytic solutions of Eqs. (15 1 and 
stage of inflation, where x 1, and in the post-inflationary 
stage, where x ~ 1- In Subsection |IIIB| we use numerical 
calculations to illustrate our analytic studies. 



A. Analytic solutions 

We first study the early stage of inflation. When x ^ 1, 
it yields that uj = luq. So Eqs. ( [T5] l and ( [T6] l return to the 
equations of motion in the extended inflation 1 17J, and we then 
obtain the power-law solutions 



a{t) = as (1 + 
xit)=XB (l + ^^)'' 



(17) 
(18) 



with fit) = 
3)(6wo + 5) and 7^ 



(1 + ^t), where 



Ml 



Here, as and xb denote the 



18 I indicates that x has a 



3a2X 

initial values of a{t) and x(^)- Eq 

constant acceleration. If > |, we obtain a power-law in- 
flation, which yields a time-dependent bubble nucleation rate 
e{t). As mentioned in [171, the initial bubble nucleation rate 
cb can be small and then grows to a critical value €cr, where 
the system of bubbles will percolate at some finite time. It 
means that e will reach tcr in the post-inflationary stage. The 
constraint of the bubble-size distribution required cjq < 25 in 
the extended inflationary model |19j,i20j|, so we may require 
w(x) < 25 in this power-law inflationary stage. Moreover, 
since ljj{x) becomes large in the post-inflationary stage, we 
should restrict the e-folding number N{t) = In to be 
less than 55 at the post-inflationary epoch in order not to pro- 
duce a large-w, scale-invariant bubble spectrum. In |26|, Lid- 
die and Wands analyzed the intermediate inflationary model, 
which has ti;($) = cjq + i^m[-S^-)™', and obtained a con- 

straint on cjq and m. They concluded that the choice of w(x) 
must exhibit a prolonged flat region and only increase rapidly 
once X approaches to 1. It corresponds to choose a large m. In 
this extended inflation model with torsion, we find that „.?^ ^ 
changes very rapidly when x approaches to 1 and only be- 
comes significant when x is extremely close to 1. So the con- 
straint of the bubble-size distribution can be achieved in this 
inflationary model by requiring wq < 20. A more detailed 
study of the bubble spectrum in this extended inflation model 
with torsion may lower the upper bound of . 

In the post-inflationary stage, we try to find an attractor so- 
lution, which asymptotically approaches to 1. More precisely, 
the solution satisfies lim^^oo X = 1 and limt^.oo = 0, 



Vn > 1. When x '■ 
and 



l,the 



-^00 X = i anc 
Eqs. ([l5|-([l6| 



become _ff ^ 



X + iHx = - 



2(1 -X) 



4M|(1 - x) 
3a2 



which yields an approximate analytic solution 



a{t) cx sinh/3t, x(^) ~ tanh^ j3t, 



3x02 



(19) 



(20) 



with f{t) 



-f3 tanh /3t, where /3 



/3a2 ' 



We see that x 



asymptotically approaches to 1, and a{t) becomes nearly ex- 
ponential expansion in this post-inflationary stage. It means 
that oj{x) will grows to a large value at the end of inflation. 
After the end of inflation, the Universe may be thermalized by 
bubble collisions and returns to radiation domination, so the 
false vacuum energy in Eq. ( [T9| should be replaced by p — 3p, 
which is zero at radiation domination. From Eq. ( [T9] l, one can 
actually see that either in the radiation- or matter-dominated 
era, the matter fields {p — 3p) do not affect the evolution of x 
since they all multiply a very small value 1 — x- It turns out 
that X will continuously approach to 1 at the radiation- and 
matter-domainated epochs and it naturally gives an extremely 
large value of uj at the present matter-dominated epoch, which 
satisfies the solar system observations. In Subsection |III B| we 
apply numerical calculations to study the evolution of x and a 
during the inflation. 
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B. Numerical demonstration 



In this subsection, we use a numerical method to demonstrate our analytic study. In the numerical calculation, we normalize 
a2 — 1 and choose Mp = 10~^ Il32l . Moreover, we set the BD parameter cjo = 16, and the initial values are chosen as as = 1, 
Xb = 10"^ and xb = 2 x 10~^. In Fig. [l] one can clearly see that xi^) is proportional to at early stage of inflation, 
which agrees with our analytic solution, and will then pass a critical point (T w 2700), where its acceleration x{T) becomes 
deceleration. In the post-inflationary stage (T > 3500), both the velocity x{T) and acceleration x{T) approach to zero, so it 
yields that x asymptotically approaches to 1, which also agrees with our analytic result. In the hiu — T diagram, we observe 
that ui is nearly ojq and grows rapidly to a large value at the post-inflationary stage. 



FIG. 1: Evolution of a{T), x{T) and oj during the inflation, where T is a dimensionless time-parameter normalized by T = 10~^Mf t. (1) 
The top left-panel indicates Ina — T diagram and the top right-panel denotes x{T) — T diagram. (2) The bottom left-panel indicates the 
evolution of the velocity of x and the bottom right-panel shows the evolution of In uj with respect to T. 
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IV. CONCLUSION AND DISCUSSION 

We study the BD theory with torsion and obtain a general 
torsion solution, which is completely determined by the di- 
mensionless BD scalar field x- We further discover that tor- 
sion fields will contribute to the BD parameter wq to form an 
effective BD parameter uj{x) — i-^o + 2(T^x)' extended 
inflation model, the constraint of bubble-size distribution re- 
quires u! < 25. However, the current solar system observa- 
tion requires cj > 40000. This apparent conflict is called the 
cj-problem. We study the evolution of x(t) and a{t) during 
the inflation, i.e false vacuum energy domination, and show 
that Li>{x) is approximate to ljq during the inflation and will 
rapidly transit to a large value in the post-inflationary stage. 
Moreover, since x continuously approaches to 1 during the 
radiation- and matter-dominated epochs, uj{x) will become 
extremely large at present time, which naturally explains the 
solar system observations. 

In this work, we solve the oj-problem by generalizing the 
BD theory to RC space-times. The next question is to un- 
derstand whether this extended inflation model with torsion 
will satisfy CMB anisotropic observations. In particular, the 
superhorizon-scale anisotropic spectrum {I < 100 modes) of 



CMB, which contains the information of primordial quantum 
fluctuations, has been used to test and constraint the inflation 
models. From WMAP 7-year data |33|, the best-fit cosmolog- 
ical parameters give a spectral index of density perturbation 
n w 0.96, which is nearly scale-invaiiant. Refs. 1341135 1 indi- 
cate that the extended inflation with w < 25 yields n < 0.85, 
which is inconsistent with the WMAP 7-year data. Further- 
more, the scalar-tensor ratio in the extended inflation with 
a; < 25 is too large to satisfy the same data |35|. In order 
to answer whether this extended inflation model with torsion 
satisfies the CMB anisotropic spectrum, further study on cos- 
mological perturbation in RC spacetimes will be followed. 
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